Acceleration of regional, racial and gender disparities in drug overdoses across the United States driven by COVID-19 and fentanyl use (preprint)

Objectives. To examine trends in drug overdose deaths by gender, race and geography in the United States during the period 2013-2020. Methods. We used the final National Vital Statistics System multiple cause-of-death mortality files to extract crude rates by gender and race and to calculate the male-to-female ratios of crude rates between 2013 and 2020. We established 2013-2019 temporal trends for four major drug types: psychostimulants with addiction potential (T43.6, such as methamphetamines); heroin (T40.1); natural and semi-synthetic opioids (T40.2, such as those contained in prescription pain-killers); synthetic opioids (T40.4, such as fentanyl and its derivatives) through a quadratic regression and determined whether changes in the pandemic year 2020 were statistically significant. We also identified states, race and gender categories most impacted by drug overdoses. Results: Nationwide, the year 2020 saw statistically significant increases in overdoses for all drug categories except heroin, surpassing predictions based on 2013-2019 trends. Crude rates for Blacks of both genders surpassed those for Whites for fentanyl and psychostimulants in 2018 creating a gap that widened through 2020. In some regions mortality among Whites decreased while overdose rates for Blacks kept rising. The largest 2020 mortality statistic is for Black males in the District of Columbia, with a record 134 overdoses per 100,000 due to fentanyl, 9.4 times more than the fatality rate among White males. Male overdose crude rates in 2020 remain larger than those of females for all drug categories except in Idaho, Utah and Arkansas where crude rates of overdoses by natural and semisynthetic opioids for females exceeded those of males. Public Health Implications: Drug prevention, mitigation and no-harm strategies should include racial, geographical and gender-specific efforts, to better identify and serve at risk groups.

A statistical model of COVID-19 testing in populations: effects of sampling bias and testing errors (preprint)

We develop a statistical model for the testing of disease prevalence in a population. The model assumes a binary test result, positive or negative, but allows for biases in sample selection and both type I (false positive) and type II (false negative) testing errors. Our model also incorporates multiple test types and is able to distinguish between retesting and exclusion after testing. Our quantitative framework allows us to directly interpret testing results as a function of errors and biases. By applying our testing model to COVID-19 testing data and actual case data from specific jurisdictions, we are able to estimate and provide uncertainty quantification of indices that are crucial in a pandemic, such as disease prevalence and fatality ratios.

A statistical model of COVID-19 testing in populations: effects of sampling bias and testing errors (preprint)

We develop a statistical model for the testing of disease prevalence in a population. The model assumes a binary test result, positive or negative, but allows for biases in sample selection and both type I (false positive) and type II (false negative) testing errors. Our model also incorporates multiple test types and is able to distinguish between retesting and exclusion after testing. Our quantitative framework allows us to directly interpret testing results as a function of errors and biases. By applying our testing model to COVID-19 testing data and actual case data from specific jurisdictions, we are able to estimate and provide uncertainty quantification of indices that are crucial in a pandemic, such as disease prevalence and fatality ratios.

Using excess deaths and testing statistics to improve estimates of COVID-19 mortalities (preprint)

Factors such as non-uniform definitions of mortality, uncertainty in disease prevalence, and biased sampling complicate the quantification of fatality during an epidemic. Regardless of the employed fatality measure, the infected population and the number of infection-caused deaths need to be consistently estimated for comparing mortality across regions. We combine historical and current mortality data, a statistical testing model, and an SIR epidemic model, to improve estimation of mortality. We find that the average excess death across the entire US is 13% higher than the number of reported COVID-19 deaths. In some areas, such as New York City, the number of weekly deaths is about eight times higher than in previous years. Other countries such as Peru, Ecuador, Mexico, and Spain exhibit excess deaths significantly higher than their reported COVID-19 deaths. Conversely, we find negligible or negative excess deaths for part and all of 2020 for Denmark, Germany, and Norway.

Using excess deaths and testing statistics to improve estimates of COVID-19 mortalities (preprint)

Factors such as non-uniform definitions of mortality, uncertainty in disease prevalence, and biased sampling complicate the quantification of fatality during an epidemic. Regardless of the employed fatality measure, the infected population and the number of infection-caused deaths need to be consistently estimated for comparing mortality across regions. We combine historical and current mortality data, a statistical testing model, and an SIR epidemic model, to improve estimation of mortality. We find that the average excess death across the entire US is 13$\%$ higher than the number of reported COVID-19 deaths. In some areas, such as New York City, the number of weekly deaths is about eight times higher than in previous years. Other countries such as Peru, Ecuador, Mexico, and Spain exhibit excess deaths significantly higher than their reported COVID-19 deaths. Conversely, we find negligible or negative excess deaths for part and all of 2020 for Denmark, Germany, and Norway.

Why estimating population-based case fatality rates during epidemics may be misleading (preprint)

Different ways of calculating mortality ratios during epidemics can yield widely different results, particularly during the COVID-19 pandemic. We formulate both a survival probability model and an associated infection duration-dependent SIR model to define individual- and population-based estimates of dynamic mortality ratios. The key parameters that affect the dynamics of the different mortality estimates are the incubation period and the length of time individuals were infected before confirmation of infection. We stress that none of these ratios are accurately represented by the often misinterpreted case fatality ratio (CFR), the number of deaths to date divided by the total number of infected cases to date. Using available data on the recent SARS-CoV-2 outbreaks and simple assumptions, we estimate and compare the different dynamic mortality ratios and highlight their differences. Informed by our modeling, we propose a more systematic method to determine mortality ratios during epidemic outbreaks and discuss sensitivity to confounding effects and errors in the data.

Why case fatality ratios can be misleading: individual- and population-based mortality estimates and factors influencing them (preprint)

Different ways of calculating mortality ratios during epidemics have yielded very different results, particularly during the current COVID-19 pandemic. We formulate both a survival probability model and an associated infection duration-dependent SIR model to define individual- and population-based estimates of dynamic mortality ratios. The key parameters that affect the dynamics of the different mortality estimates are the incubation period and the time individuals were infected before confirmation of infection. We stress that none of these ratios are accurately represented by the often misinterpreted case fatality ratio (CFR), the number of deaths to date divided by the total number of confirmed infected cases to date. Using data on the recent SARS-CoV-2 outbreaks, we estimate and compare the different dynamic mortality ratios and highlight their differences. Informed by our modeling, we propose more systematic methods to determine mortality ratios during epidemic outbreaks and discuss sensitivity to confounding effects and uncertainties in the data.